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Abstract. The recent interest in using microorganisms for biofuels is motivation 
enough to study bioconvection and cell dispersion in tubes subject to imposed flow. 
To optimize light and nutrient uptake, many microorganisms swim in directions bi- 
ased by environmental cues (e.g. phototaxis in algae and chemotaxis in bacteria). Such 
taxes inevitably lead to accumulations of cells, which, as many microorganisms have 
a density different to the fluid, can induce hydrodynamic instabilites. The large-scale 
fluid flow and spectacular patterns that arise are termed bioconvection. However, the 
extent to which bioconvection is affected or suppressed by an imposed fluid flow, and 
how bioconvection influences the mean flow profile and cell transport are open ques- 
tions. This experimental study is the first to address these issues by quantifying the 
patterns due to suspensions of the gravitactic and gyrotactic green bifiagellate alga 
Chlamydomonas in horizontal tubes subject to an imposed flow. With no flow, the 
dependence of the dominant pattern wavelength at pattern onset on cell concentration 
is established for three different tube diameters. For small imposed flows, the vertical 
plumes of cells are observed merely to bow in the direction of flow. For sufficiently high 
flow rates, the plumes progressively fragment into piecewise linear diagonal plumes, un- 
expectedly inclined at constant angles and translating at fixed speeds. The pattern 
wavelength generally grows with flow rate, with transitions at critical rates that de- 
pend on concentration. Even at high imposed flow rates, bioconvection is not wholly 
suppressed and perturbs the flow field. 

Keywords: swimming micro-organisms, algae, Chlamydomonas, bioconvection, gyro- 
taxis, photobioreactors, pipe flow. 
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1. Introduction 

Most populations of motile microorganisms bias their motion in response to 
environmental cues to seek better conditions (e.g. bacterial chemotaxis, Berg 2004). To 
find water surfaces, phototrophic green algae, such as Chlamydomonas and Dunaliella, 
are gravitactic (they swim upwards on average). Typically the upward bias is due to a 
combination of the cell's centre-of-mass being displaced from its centre-of-buoyancy and 
a torque due to sedimentation of the asymmetric body-flagella ensemble. Additionally, 
such cells are often gyrotactic (on average they swim towards downwelling regions due 
to a balance of viscous and gravitational torques; Kessler 1985) and phototactic (biased 
motion relative to the illumination; Foster & Smyth 1980). The taxes invariably lead to 
accumulations of the negatively buoyant cells, inducing hydrodynamic instabilities that 
result in spatially localized structures and patterns, termed bioconvection (see Pedley 
& Kessler 1992; Hill & Pedley 2005). In particular, overturning instabilities occur when 
cells accumulate at the upper surface, forming a dense, unstable layer (analogous to 
Rayleigh-Bernard convection; Chandrasekar 1961), yielding patterns of spots, stripes or 
more complicated structures. Away from the upper boundary, gyrotactic instabilities 
lead to long thin plume structures, where cells swim towards downwelling regions and 
their added mass amplifies the downwelling. 

Recently, Durham et al. (2009) showed that transient gyrotactic trapping of 
motile phytoplankton in differentially sheared fluids can generate thin, cell-rich 
layers in the ocean, with important consequences for the ecology of toxic algae. 
Other investigations (e.g. Bearon & Griinbaum 2006; Grunbaum 2009) demonstrate 
the ecological importance of spatially localized structures (patchiness, due to cell 
swimming). In biotechnology applications, the properties of swimming cell suspensions 
have been minimally exploited or outright ignored. Recently, there has been renewed 
interest in the production of biofuel from microorganisms (especially biodiesel or 
hydrogen; Melis & Happe 2001; Chisti 2007). For these alternative fuels to become 
commercially competitive, existing bioreactor designs need to be refined. In current 
closed bioreactors, algae are cultured in laminar or turbulent flows through arrangements 
of horizontal, slanted and/or vertical tubes, bubbled for gas exchange, and concentrated 
by flocculation, filtration or centrifugation, as if they were chemicals or inert colloids 
(Grima et al. 2001, 2003; Garcia- Gonzalez et al 2005; Chisti 2007). However, candidate 
species for biofuel production that swim (e.g. Dunaliella) do not behave like chemicals 
or inert colloids: their peculiar collective motions and consequent transport properties 
need to be taken into account, and indeed exploited. Whether the flow in tubes is 
laminar or turbulent, gyrotactic swimming cells organize in patterns that alter the flow. 

Bees & Croze (2010) analysed the dispersion of swimming cells in a flow in a 
vertical tube in the laminar regime. They predicted the effective drift and diffusion of 
gyrotactic algae in plumes; depending on the algal cell properties, algal dispersion can 
be quantitatively very different from that of classical Taylor-Aris dispersion for inert 
chemicals. Bees & Croze (2010) also obtained steady state plume solutions, which show 
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Figure 1. Bioconvection pattern of C. augustae in a tube (tube A; concentration, 
c « 1.4 X 10 6 cells/ml) 8 minutes after the onset of instability, as imaged: (a) from 
the top, focusing close to the top of the tube; (b) from the side, focusing on the center 
of the tube. The observed patterns are similar to those originally photographed by 
Wager (1911) for E. viridis. The images are 4.2 by 1.3 cm. 

that the flow in the presence of cells deviates from the standard Poiseuille profile. 

It is well known that the transition to turbulence in a tube is strongly sensitive to the 
initial laminar state (Willis et al 2008). It is interesting to note that the perturbation by 
the presence of swimming cells should inevitably change the onset of the transition. Cells 
in turbulent flows in pipes have not been analysed, but Lewis (2003) recently showed 
that gyrotactic algae in a homogenous and isotropic turbulent flow field retain their 
bias. Turbulence only changes the effective value of the diffusivity of cell orientation. 

In this paper, we present an experimental exploration of bioconvection in horizontal 
tubes under conditions of no imposed flow and in the presence of weak shear. We aim 
to investigate the transport of algae in a minimally perturbed system with bioreactor 
geometry, in the same spirit as Bees Sz Croze (2010), and with a view to understanding 
the coupled cell and fluid dynamics in bioreactors. This is the first study of bioconvection 
under shear, and the first quantitative study of bioconvection in horizontal tubes. The 
only other experimental study of bioconvection in horizontal tubes is the pioneering 
work by Wager (1911). Wager placed suspensions of the alga Euglena viridis in 25 cm 
long tubes of diameter 0.7 cm, and described and imaged the resulting 'aggregations' 
of cells (see figured]): an instability was visible after 30 s and the patterns evolved to a 
stable state over tens of minutes; cell concentrations were not established, but patterns 
were sharper and formed quicker for more concentrated samples. Wager repeated the 
experiments with Chalmydomonas cells (20 cm long tubes of diameter 0.8 cm), obtaining 
similar structures (but requiring larger concentrations). 

There are few quantitative experimental studies of bioconvection and these have 
been limited to shallow layers with a free upper surface. Bees & Hill (1997) 
analysed bioconvection patterns of Chlamydomonas nivalis (in this paper denominated 
Chlamydomonas augustae, see methods) in 5 cm diameter Petri dishes. From Fast 
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Fourier Transforms (FFT) of processed pattern images they obtained the concentration, 
depth and time dependence of the initial and long-term pattern wavelengths. Both 
wavelengths were found to fall off with concentration. The initial wavelength also 
increased with depth (see Czirok et al 2000), while the final wavelength did not strongly 
depend on it (due to gyrotactic focusing). Using a similar approach, the wavelengths 
of the bioconvection patterns due to the aerotactic bacterium B. subtilis have been 
quantified (Janosi et al 1998; Czirok et al 2000). The initial pattern wavelength 
depended on cell concentration but not depth. The differences between species were 
ascribed to the particular taxes in operation: aerotactic bacteria are in thin layers where 
oxygen gradients exist, whilst gravity affects gyrotactic cells throughout the layer. 

Theoretical studies of bioconvection are more numerous. For green algae, 
continuum models of the coupled fluid flow and cell concentration have evolved 
significantly from deterministic models of pure gravitaxis (Levandowsky et al 1975, 
Childress et al 1975) to stochastic models of gyrotactic cells (Pedley & Kessler 1990). 
Predictions for the wavelengths at the onset of bioconvection from linear stability 
analyses give values of the same order of magnitude as experiments: slightly larger for 
uniform suspensions (Pedley & Kessler 1990) and smaller for layers of finite depth (Bees 
& Hill 1998). The discrepancies between theory and experiment have been attributed 
to the uncertainty with which key experimental parameters are known, nonlinear effects 
and to the non-attainment of equilibrium before instability arises [idem). Nonlinear 
effects have been explored by Bees & Hill (1999) who calculated static and travelling 
plume solutions for bioconvection in deep layers. We refer the reader to the reviews by 
Pedley & Kessler (1992) and Hill & Pedley (2005) for further details. 

In the next section, we describe the cultures, experimental set-up and data 
acquisition methods used in this work. In section [3] we present results quantifying 
bioconvection patterns first in the absence and then in the presence of imposed flow. 
Finally, in section H] we discuss these results and draw conclusions from our observations. 

2. Materials and methods 

Cell culture and concentration 

Batch cultures of C. augustae (CCAP 11/51B; Culture Collection of Algae and Protozoa) 
were grown statically in conical flasks to exponential phase on 3N-BBM (nitrogen- 
enriched Bold's basal medium; Schlosser 1997) at a stable ambient temperature of 
T = 24 ± 1°C. The up-swimming cells accumulate on top of sterile cotton wool placed 
in the neck of a full flask and can be harvested by gentle suction with a Pasteur pipette 
after 2 or 3 days, and diluted to provide a range of concentrations. Concentration 
was measured with a spectrophotometer (WPA CO7500) from readings of A 590 , the 
absorbance at 590 nm, calibrated with a haemocytometer. Note that the green alga 
known as C. nivalis (CCAP 11/51B, see Kessler et al (1992)), the model organism for 
bioconvection studies, was recently shown to have been mistakenly identified (Proschold 
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Table 1. Table of diameters and lengths of different cylindrical tubes used in the 
experiments. Tubes are made of pyrex glass except tube C which is plastic. 



tube 


inner/outer diameter (cm) 


length (cm) 


A 


1.10 / 1.35 


30.0 


B 


0.79 / 1.03 


28.5 


C 


0.59 / 0.81 


33.0 




Figure 2. Schematic side-view of the experimental setup. The center of the tube is 
encased with a glycerol-hllcd flat-faced perspex bottle to minimize optical distortion, 
clamped where indicated by the shaded boxes and illuminated by a red LED light 
source fixed behind the tube. A camera (not shown) mounted in front of the tube 
captures images from the side. Suspensions are loaded via a plastic tube affixed to the 
left of the tube. Flow (in the direction of the arrow) is imposed using a syringe pump 
on the right of the tube, where we also see a tap to facilitate air removal when loading. 

et al 2001). It has been renamed C. augustae, and we use this denomination herein. 
Set-up, mixing, illumination and flow 

Three tubes of different diameters (pyrex glass or plastic; table CP were employed (see 
figure |2]) . The middle of each tube was encased with a flat-faced perspex bottle filled 
with water or glycerol to minimize optical distortion. Results using glycerol, where 
distortion is small but non-negligible at the tube edges, are systematically corrected 
(see below). Tubes were clamped about the central casing, placed on an isolated work 
bench to minimize vibrations, and levelled before introducing the suspension. 

Tubing was attached to the ends of the cylindrical tubes to aid mixing and to 
drive flows (figure [2]). Suspensions were mixed in a plastic bottle using a slow magnetic 
stirrer before being poured into the tube and imaged. The pouring (and any inadvertent 
entrainment of bubbles) generated secondary flows, and appeared to aid mixing. Other 
mixing protocols gave less satisfactory results. The proportion of swimming cells was 
reduced by repeated mixing due to deflagellation and adhesion of cells to the boundaries. 
To prevent the latter, we tried soaking the tubes in bovine serum albumen (BSA), 
which limits adhesion in microchannels (Weibel et al 2005), but patterns were slightly 
suppressed. The decrease of the concentration by surface fouling was not significant 
on the experimental timescale (10 minutes). It is conceivable that in the longer lasting 
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experiments with flow, the flow itself helped prevent adhesion. 

Tubes were illuminated with diffuse red light of wavelength 660 nm (BL1960 LED 
light bank; Advanced Illumination, Rochester, USA; at less than 100 lux), to which 
green algae do not respond phototactically (Nultsch et al 1971; Foster & Smyth 1980). 

Flow in the tubes was generated using a Graseby 3500 syringe pump (Graseby 
Medical Ltd., Watford, UK) attached at one end of the experimental setup (figure [2]) 
to induce a pressure gradient that results in a weak Poiseuille flow if cells are absent. 
Experiments were performed with flow rates in the range — 35 ml/h, as described in 
detail the next section, at a stable ambient temperature of T = 24 ± 1°C. 

Data acquisition, processing and analysis 

Timelapse sequences of projected bioconvection patterns were captured from the side 
and from above using either a computer controlled CL-1014 B/W CCD camera (Camtek) 
or EOS 350D SLR digital camera (Canon) for higher resolution results. In the absence 
of flow, up to 201 images were recorded at intervals of 2 to 5 s, whilst for the flow 
experiments 49 images were recorded a second apart. Images were acquired (768 x 576 
pixels for CL-1014; 3456 x 2304 pixels for EOS 350D) and processed using IDL (RSI, 
Boulder, USA). The first image of a sequence, if featureless, was subtracted to remove 
static noise; when features appeared in the first image (flow results), the stack median 
was subtracted instead. Each cropped image, h(x,y), was windowed before applying 
the FFT to give H(k x ,k y ) (see Bees & Hill 1997), yielding the biased spectrum 

/(*)= i#(^A)i 2 , (i) 

d=k,\i>\<w/4: 

where d(k x , k y ) = ^Jk^ + ky and if)(k x , k y ) = axctaii(k y /k x ) is the angle to the horizontal. 
In (pQ), Fourier modes (k x ,k y ) with \ifj\ > it /4 were filtered to eliminate contributions 
from the boundaries of the tube and cross-sectional structure. The double logarithmic fit 
proposed by Czirok et al (2000) was employed: ln[/(/c)] = a\ ln(k) — ln(k ) \ — (3 ln(&) + c, 
where a, (3, c and ko (the dominant wavenumber) are fitting parameters. Wavelengths 
are given by A = Iw/k , where Iw is the image width. The initial wavelength, Aj, 
was determined from the fitted Fourier spectra together with direct observation of 
image sequences and contour plots (figure [3}) as a guide when bimodal distributions 
(occasionally) occurred. The time for the establishment of the initial instability was 
measured as the time between pouring in a mixed suspension and measurement of the 
initial wavelength. The final pattern wavelength, A/, was calculated as an average over 
the last 50 images of each sequence; the standard deviation of this sample was taken as 
an estimate of uncertainty in the value of A/. When more than one value was obtained 
for a given concentration, the weighted average and standard error in the mean of the 
results were reported. 

With imposed flow, images were acquired as follows. The cell suspension was mixed 
and slowly poured in as before. Then, after waiting 3 — 5 minutes for the pattern to 
form, 49 images were recorded. The flow rate was then increased to 2 ml/h, and after 
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waiting 2 minutes, another 49 frames were recorded. Proceeding similarly, the flow 
rate was increased in steps of 2 ml/h up to 10 ml/h, and then in steps of 5 ml/h 
up to a maximum of 35 ml/h. (For some experiments only 5 ml/h steps were used). 
The dominant wavelength was obtained by averaging over the last 20 frames of each 
sequence (the standard deviation was taken as an estimate of the variation). To quantify 
pattern distortion and plume dynamics, the mean angle the plumes make to the vertical 
towards the top and bottom of the tube, (9), and the mean plume drift speed, (V), 
were measured. Plume angle measurements were obtained using a straight line selection 
tool in ImageJ, NIH, for the approximately straight sections of the plumes. Straight 
plume fragments were measured directly in the upper or lower regions of the tube or, 
when curved plumes spanned the whole tube, measurements were taken at the top and 
bottom boundaries. Angles were averaged across plumes from three clear images in each 
sequence (the standard deviation was used as a measure of uncertainty). For the plume 
drift speeds, each plume was tracked in a sequence and its displacement, L p , recorded 
over a total time T p , giving the mean plume speed, V = L p /T p (the portion of the plume 
that does not change shape significantly was tracked; the plumes were not observed to 
accelerate). The mean of the speeds over different plumes provides (V). 

For results with the glycerol filled jacket, the image distortion due to refraction at 
the water-glycerol interface was corrected: if z' is the measured vertical position of a 
pixel, z is the actual vertical position (z — at the tube center) and v = 1.10 is the 
ratio of the index of refraction of pyrex/glycerol to that of water, then approximately 
\z\ ~ v\z'\. Distortion only operates vertically, and marginally affects the measurement 
of angles as described above (corrected via tan((#)) tan( (#'))/ v\ less than a 3° change). 

3. Results 

Here we report the results of our observations and quantitative analysis of bioconvection 
in horizontal tubes, first in the absence of flow and then for sheared bioconvection. 

3.1. No flow 

In the absence of flow, the pattern development was recorded after a cell suspension 
was mixed and poured into each tube. Experiments were carried out with tubes A, 
B and tube C (inner diameters 1.10,0.79 and 0.59 cm respectively), varying the cell 
concentration, c, in the range c = 0.5 x 10 6 — 4 x 10 6 cells/ml. The sequence in figure 
|3]^a-h) shows the evolution of a typical bioconvection pattern in a tube as viewed from 
the side (tube A, c = 1.84 x 10 6 cells/ml). We observe that a well mixed suspension 
(homogeneous in the first frames) quickly becomes unstable, breaking up into a beautiful 
striated alternation of thin dark plumes (high cell concentration) and white bands (low 
cell concentration). Similarly to bioconvection patterns in other geometries, the pattern 
is seen to sharpen in time, with plumes becoming spaced closer together. As one 
would expect for cylindrical tubes of uniform diameter, statistically similar patterns 
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15 20 
Wavenumber 



Figure 3. Image sequence (a-h) and corresponding time- wavenumber contour plot 
(i) of the Fourier spectrum density of a representative bioconvection pattern (tube A, 
c = 1.4 x 10 6 cells/ml). In the sequence the processed images are shown at times 20, 
35, 45, 60, 75, 160, 260, 435 s (a-h). The sequence and contour chart the evolution of 
the pattern from the onset to the final steady state. 



are observed along the length of the tube. 

The above observations are usefully quantified by analysis of the dominant Fourier 
modes extracted from the images (see methods). Figure El shows a contour plot of the 
power spectrum of the image sequence in figure [3]^ a-h), demonstrating the evolution of 
pattern wavenumbers with time. The Fourier spectrum quantifies our observations: an 
initially homogenous suspension quickly suffers an initial instability, with wave number 
ki (the first 'peak' in figure [3}); the pattern then becomes unstable to other modes, 
evolving to a structure with an average final wave number kf. The complex dynamics 
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Figure 4. Initial (a) and final (b) dominant wavelengths, \i and A/ respectively, as a 
function of concentration for tubes A, B and tube C (top three panels, as indicated). 
Straight line fits to the data are shown (for tube C in (a) the fit was for c > 1.5 x 10 6 
cells/ml, see text). The last panel shows the fits without the data to allow a comparison 
of magnitudes. The error bars shown are standard deviations over a sample of 50 fits 
in a sequence (b) or the standard error in the mean over repeated experiments (a, b). 



with which modes evolve from the initial instability is also clear in figure [3}. From 
figure HI it is clear that the initial and final dominant pattern wavelengths, \ and A f 
respectively, are consistent with a linear drop with concentration. A possible exception 
is the narrow tube C, where the large rise in Aj at very low concentrations might be 
evidence for stronger dependence. Thus, excluding points below c = 1.5 x 10 6 cells/ml 
in tube C for Aj, the data can be described by the linear form A(c) = Ao + crc, where 
Ao (intercept) and a (gradient) are the fit parameters. Performing a least squares fit 
for the initial wavelength: (of = -0.085 ± 0.030 cm ml/cells, A^ = 0.84 ± 0.08 cm) 
for tube A; (of = -0.093 ± 0.023 cm ml/cells, A^ = 0.64 ± 0.06 cm) for tube B; and 
(of = -0.075 ± 0.028 cm ml/cells, A^ = 0.70 ± 0.08 cm) for tube C. For the final 
wavelength we find similar values, (of = —0.076 ± 0.049 cm ml/cells, A4 = 0.65 ± 0.13 
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Figure 5. Variation of time to establish the initial instability for different cell 
concentrations for tubes A, B and C (top three panels, as indicated), and comparison 
of linear fits to this variation (bottom panel) . 



cm) for tube A; (exj = -0.046 ± 0.013 cm ml/cells, Af = 0.45 ± 0.03 cm) for tube B; 
and (af = -0.088 ± 0.004 cm ml/cells, A$? = 0.57 ± 0.01 cm) for tube C. Within the 
uncertainties of our measurements, it is hard to ascertain any dependence of wavelength 
on tube diameter, but it would be reasonable to infer that both the initial and final 
wavelengths are in general larger for tube A than for tube B or C. 

Aside from the initial wavelength, another useful measure of the onset of 
bioconvection is the time Tj taken to establish the initial instability (figure [5]). Like 
Aj, Tj falls linearly with concentration, but without any obvious sudden rise at low c for 
tube C. Using a least squares fit with a function of the form Tj(c) = r ti + jiC, (with r j 
and ji fit parameters, as above), we find: (7^ = —6.0 ±2.1 s ml/cells, = 43 ±6 s) for 
tube A; (7? = -5.6 ± 1.9 s ml/cells, r£ = 38 ± 5 s) for tube B; and (7? = -6.4 ± 2.7 s 
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Figure 6. Bioconvection patterns in tube B subject to imposed flow in the range 
Up = to 16.6 cm/s (F = to 35 ml/h) for cell concentrations (a) c = 1.26 x 106 
cells/ml, (b) c = 2.22 x 106 cells/ml and (c) c = 5.18 x 106 cells/ml. (As indicated, 
the second frame in (a) shows a frame for 3.4 cm/s rather than 2.8 cm/s). (d) For 
Up = 13.8 cm/s and concentration as in (b), we show three 5 s interval snapshots of the 
dynamics of a bowed plume (indicated by arrows). The diagram to the right displays 
tracings of the plume (scale bar is 0.2 cm). 



ml/cells, = 41 ± 7 s) for tube C. The dependence on tube radius again is not strong. 
3.2. Effect of flow 

When flow was imposed using the syringe pump, patterns could be observed to be 
progressively distorted. Figure |6] shows the patterns in tube B for (a) c = 1.26 x 10 6 
cells/ml, (b) c = 2.22 x 10 6 cells/ml and (c) c = 5.18 x 10 6 cells/ml as the mean flow 
speed Up is increased in steps of 2.8 cm/s from no flow to a maximum of 16.6 cm/s 
(the flow is from right to left). The qualitative behaviour is generally similar for the 
different concentrations studied. In the absence of flow, the bioconvection plumes are 
parallel (on average) and slightly denser at the bottom. For small to moderate flow 
rates, most plumes still span the tube, but are seen to be curved so that the top and 
bottom of a plume are tilted at a characteristic angle (on average) to the vertical. The 
bowed plume shape is not symmetric for low concentrations, with plumes originating 
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at the top stretching beyond the tube midpoint, before bending back at the bottom. 
Observation of the flowed pattern sequences reveals that plumes translate horizontally 
with the flow and that often high concentration pulses (blips) travel down a plume. As 
the flow rate is increased further, the pattern gets distorted more and more, with plumes 
tilting further away from the vertical at the top and bottom. The pattern appearance for 
higher flow rates is less orderly; many plumes break up, with less plumes spanning the 
tube. However, at high flow rates, the plumes becomes more symmetric: linear bottom 
standing plumes stretch to the tube midpoint. Remarkably, the plumes maintain a 
constant average angle at the top and bottom of the tube, which increases with flow 
rate, as plumes are tilted by the flow (see below). The reasons for this are nontrivial, see 
discussion. It should be remarked that the plume behaviour reported above is statistical: 
sequences show that individual plumes are not stable, meandering and either focusing or 
dispersing. For the highest flow rates all plumes fragment into very elongated structures 
making a large angle to the vertical and top and bottom plumes become fully symmetric 
about the tube midplane. As stated above, the phenomenology described is qualitatively 
similar for all concentrations studied, although there are quantitative differences. 

To quantify the inclination of plumes, we measured the average angle, (9), a plume 
makes to the vertical when measured from the top and bottom of the tube (see methods). 
For top plumes, 9 is defined as positive for tilt away from the vertical in the direction 
of the flow, while for bottom plumes it is negative. We see from figure [7b, that plume 
tilt does not seem to vary sensitively on concentration. The apparent independence of 
(9) on concentration can be made clearer by plotting all concentrations on the same 
graph, as in figure [H^i, where the change in (9) with flow rate F has been fitted with a 
linear trend constrained to go through the origin (since we expect (9) = for F = 
ml/h). The gradients obtained are very close in absolute value: 3.65 ±0.11 degrees 
s/cm for the top plumes and —3.91 ± 0.13 degrees s/cm for the bottom plumes. The 
average drift speed, (V), of the plumes was also measured. As shown in [8b, where all 
concentrations are plotted on the same graph, (V) appears to increase linearly with flow 
rate in a manner which appears generally independent of concentration. A constrained 
linear fit to the data, yields a value 1.30 ± 0.04 for the change in average plume speed 
with mean flow speed. Also shown in [8b for comparison is the linear prediction for the 
maximum speed of the Poiseuille flow, U po i = U = 2F/na 2 and for the equivalent plug 
flow U p = U/2 (i.e. the mean Poiseuille speed), giving gradients of 2 and 1, respectively. 
The data appears to fall between these two limits but closer to plug flow. We discuss 
the significance of this below. 

A plot of the dominant pattern wavelength A as a function of flow rate F (figure 
UJl) shows a clear concentration dependence of the flowed pattern dynamics. For all the 
concentrations studied, as the flow rate is increased, A initially grows with flow rate. 
However, there are critical concentrations around which this growth is interrupted by a 
sudden rise and fall of A, which quantitatively demarcates a dynamical transition in the 
pattern. These transitions represent a statistical change in the dynamics which is not 
immediately obvious from the frame sequences, however it is clear that these statistical 



Sheared bioconvection in a tube 



13 




Figure 7. (a) Final pattern wavelength and (b) average top (positive) and bottom 
(negative) plume angles to the vertical as a function of flow rate F for concentrations in 
the range c = 0.59 x 10 6 — 5.18 x 10 6 cells/ml, as indicated. The peaks in the wavelength 
around certain critical flow rates could be the signature of dynamic transitions in plume 
arrangement. The error bars shown are standard deviations. 
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Figure 8. Plot of (a) the average angles and (b) average plume drift speed as functions 
of the mean flow speed for all concentrations. Neither of these measures appears to 
depend strongly on concentration. Linear fits to data for all concentrations are also 
shown, and in the case of (b) the plume drift speed, the fit (solid line) is compared 
with the prediction for the mean (dot-dashed) and maximum (dashed) flow speeds for 
Poiseuille flow. 



changes in wavelength are associated with rearrangements of the plumes as the flow rate 
is increased. 

4. Discussion 

We have presented the first quantitative analysis of bioconvection patterns in horizontal 
tubes in the presence and absence of flow. A well mixed suspension of the green alga 
C. augustae becomes unstable to cell concentration fluctuations leading to a system 
of vertical plumes along the length of the tube. With no flow imposed, the pattern 
sharpens with time, and the plume spacing decreases, consistent with observations by 
Wager (1911) of Euglena viridis and Chlamydomonas in similar tubes. To quantify 
these observations we Fourier analysed pattern sequences imaged from the side to 
measure initial and final dominant pattern wavelengths. Both decrease slowly with 
concentration in the range 0.5 x 10 6 — 4.5 x 10 6 cells/ml (except for a steeper drop for 
the initial wavelength for low concentrations in the smallest tube). Bees & Hill (1997) 
observed a similar decrease with concentration for both the initial and final wavelength 
for bioconvection in shallow layers. They also found that the initial wavelength increased 
linearly with depth (see Czirok et al 2000), while the final wavelength appeared less 
sensitive to it. The results here hint at a similar weak dependence on tube radius for a 
fixed concentration, but this cannot be inferred with confidence. It would be interesting 
to perform experiments with a greater range of tube diameters to explore these trends. 
It should be noted that we have studied projections onto a vertical plane of what are 
clearly three-dimensional patterns (see figure 1). The cross-sectional plume and flow 
structure may have a significant impact on the above interpretations of the results. 

We have presented the first experimental study of bioconvection in the presence 
of imposed flows (with mean axial flow speeds that are comparable to bioconvective 
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circulation). We find that plumes are simply distorted for weak flows, but break as 
the flow rate increases. In this general phenomenology, the average plume drift speed 
(V) and inclination to the vertical (8), for a given flow rate, do not appear to depend 
appreciably on concentration in the range studied. However, plots of final pattern 
wavelength as a function of flow rate for different concentrations reveal that transitions 
in the pattern dynamics occur at critical flow rates that are sensitive to concentration. 
These transitions are more apparent from Fourier analysis than direct observation of 
the image sequences. 

It will be useful to discuss our results in terms of approximations of the most recent 
model of bioconvection. However, the solution of the full model is beyond the scope of 
this paper, so we shall resort to significant simplification later in this section to help 
explain our results. We briefly summarize the model here, and refer the reader to Pedley 
& Kessler (1992) for a full exposition. In a horizontal tube, the incompressible fluid flow 
velocity, u(r, t), is coupled to the cell density, n(r, t), via the Navier-Stokes equations 
with a negative buoyancy term due to the cells, whilst a cell conservation equation 
describes flow dependent cell diffusion, swimming and advection (see Bees & Croze 
2010 for a vertical tube). These flow dependent terms are obtained from a microscopic 
model of the hydrodynamics of swimming cells. 

With no imposed flow, preliminary calculations (not shown) reveal the existence of 
a steady state profile in a horizontal tube, with most cells concentrated towards the top, 
in line with bioconvection between two horizontal boundaries. Such a distribution may 
lead to an overturning instability, as in analyses of other geometries (Hill et al 1989; 
Bees & Hill 1998). Gyrotactic instabilities may also result where the cells are focused 
into plumes as a result of self-driven flow. One may question whether overturning or 
gyrotactic instabilities occur before the steady state has been approached. When cells 
are placed in the tube, fluid motions decay with a time scale, r v = a 2 / u, where a is the 
tube radius and v is the kinematic viscosity. For the largest tube (tube A) a ~ 0.55 
cm, and taking v = 10~ 2 cm 2 /s, we obtain t v = 30 s. An unstable density profile will 
be established in the time it takes a cell to swim the lengthscale a, r Q = a/V s , where 
V s is the cell swimming speed: with V s ~ 10 -2 cm/s, r = 55 s. We can compare 
these estimates with the observed time for the onset of the initial instability, Ti ~ 30 s 
(c ~ 2.5 x 10 6 ; figure [5]). We see that t v ~ r« < r ; for the largest tube, the instability 
occurs roughly when mixing effects have decayed, but a little before the steady state 
profile has been established. A similar calculation yields t v ~ r, < r D for tube B, and 
t v < Ti ~ r for tube C (the steady state has been established before pattern onset 
with negligible mixing effects). Since Ti increases for smaller c, patterns will likely be 
generated by an overturning instability for dilute concentrations. 

The observed evolution of patterns in the presence of flow results from coupling 
between the imposed suspension flux and bioconvective circulation. The flow is laminar 
for the flow rates that we study; the Reynolds number, Re = Ua/v, based on the 
maximum Poiseuille flow speed U, is 0(1). The cells are advected horizontally, but 
plumes of falling cells redistribute fluid, modifying the flow. In general, the secondary 
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flows significantly alter the mean flow profile from its cell-less Poiseuille state, and 
the results suggest that the mean profiles are sometimes closer to plug flow. For 
spherical cells, a simple balance of gravitational and viscous torques provides the angle 
to the vertical, ip, at which gyrotactic cells swim in a shear flow with vorticity u, 
such that sin^ = Bu, where B = 3.4 s is the gyrotactic reorientation time (for an 
improved model see Pedley & Kessler 1990, Bees et al. 1998; stable orientation requires 
uj < u c = 1/B = 0.3 s _1 , otherwise the cells tumble). The Poiseuille flow profile is 
u{r) = U[l — (r /a) 2 }, giving vorticity uj{r) = (2U/a 2 )r, whereas plug flow has u(r) = U/2 
and U)(r) = 0. If Poiseuille flow is assumed with mean flow rates of 0.003 to 0.02 cm 
s _1 , the largest shear is smaller than u c , with u max = a; (a) = 0.03 to 0.2 s" 1 and thus 
maximum inclinations ranging from 6° to 43°. Hence, we do not expect the formation 
of layers induced by gyrotactically trapped tumbling cells, as observed by Durham et 
al. (2009), except perhaps near the boundaries for flatter flow profiles. 

The results indicate that for small to moderate flow rates, plumes remain whole 
but are distorted by the flow, starting at the top of the tube, tilting from the vertical 
and then bending back at the bottom of the tube. This suggests that the circulation 
caused by descending plumes dominates that due to the imposed horizontal advection. 
A striking feature of the plumes is that they appear to translate horizontally at a 
fixed average speed (V), mostly preserving their shape during their lifetime. The 
progressive 'breaking' of curved plume structures into angled linear plumes as the 
flow rate is increased further is harder to interpret. Clearly stronger imposed flows 
disrupt bioconvective circulation to some degree. One hypothesis is that for small 
flows the bioconvective circulation cells adjust and translate with the mean flow profile 
speed, but when the flow becomes too strong the plumes stagger, splitting the large 
closed streamlines. The dynamical transitions highlighted from measurements of the 
dominant pattern wavelength are intriguing. From the results in figure [7] we infer that 
pattern transitions occur at higher flow rates for larger concentrations. This makes 
sense intuitively if we consider that for higher concentrations we expect a stronger cell- 
induced circulation, which will be more stable to perturbations by the imposed flow. 
The resulting picture is that patterns are affected at all nonzero flow rates, but are 
severely disrupted, with dramatic transitions in the plume dynamics, only beyond a 
concentration dependent critical flow rate. The transitions only affect bioconvective 
circulation, and so plume arrangement, which would explain why the average angle (9) 
and average plume speed (V) do not depend on concentration. 

To help gain insight into cell dynamics in flows, we shall construct a much simplified 
description. Consider the plume structure in a vertical plane at the centre of the tube, 
with cartesian coordinates (x, y) as in figure [9j We shall restrict attention to two- 
dimensional space for simplicity. At its simplest, the position of a cell is subject to 
horizontal advection by the mean flow profile, u(y), and vertical advection in a plume, 
with constant speed V. Furthermore, each cell swims at an angle to the vertical due to 
gyrotaxis with velocity V c = V s (—Buj(y), a/1 — (Buj(y)) 2 ), such that 

x(t) = u(y)-V s Bu(y) and y(t) = -V + V s ^l - B 2 u(y) 2 . (2) 
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Figure 9. Evolution of cell positions for (a) Poiseuille and (b) plug flows, according 
to the simplified description (see text). A line of cells (thin lines) released at equally 
spaced positions across the tube collapses and deforms while advected for (a), while 
it collapses vertically in (b). The stages in the evolution are 20 s apart. The upper 
envelope (thick line) is the trajectory of cells released at the top of the tube, for the 
two flow scenarios. 

Figure M presents the solution of (j2]) for different initial conditions (xq, yo), with U = 0.01 
cm/s, V = 0.0175 cm/s, V s = 0.01 cm/s, B = 3.4 s and a = 0.4 cm. The thick lines 
indicate the path of a cell released at the top of the tube (0, a), and is representative of 
a plume that is fed from the top. In Poiseuille flow (fig [9^) the path is curved due to the 
flow profile, and at the top and bottom shear induces diagonal gyrotactic reorientation 
(arrows). In plug flow (fig [9b ) the cells on average move diagonally downward but swim 
upwards. Also shown (via a series of thin lines) is the evolution of a vertical line of 
cells, initially across the tube, at time intervals of 20 s. For Poiseuille flow, the line 
translates with the horizontal component of the swimming velocity at the bottom as 
it collapses and deforms. For plug flow the line remains vertical as it collapses and 
translates. This simple description suggests that the behaviour of the plumes that 
we observe experimentally could result from a mean flow that is intermediate between 
Poiseuille and plug flows. In particular, the curved plumes are reminiscent of those for 
Poiseuille flow but the curved plumes translate with a fixed shape as though they were 
in plug flow. Furthermore, figure [8b demonstrates that plumes drift at a speed closer to 
plug flow than Poiseuille. The model demonstrates that because advection by the flow 
and gravity are not collinear, gyrotactic reorientation by the imposed flow does not lead 
per se to focusing (convergent trajectories) as in vertical tubes (Kessler 1985), although 
locally cells clearly drive the flow and will focus as a result. 

The implications of the findings for the transport of cells in tubes are of particular 
interest. In vertical tubes, cells are subject to a modified Taylor-Aris dispersion: they 
are transported within a plume parallel to the flow, with advective, swimming and 
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diffusive contributions (Bees & Croze 2010). In horizontal tubes, the results presented 
here indicate that the situation is different, with cells transported in numerous plumes 
mostly translating at fixed speeds, (V). 

To conclude, we have studied bioconvection in tubes with and without flow. 
The results demonstrate that even very weak laminar flows are sufficient to perturb 
bioconvection patterns and complex structures persist in the flow for larger flow rates, 
rather than flow simply having a mixing effect. The presence of the cells appears to 
modify the mean flow to be more like plug flow than Poiseuille. Future experimental 
investigations should measure the (mean) horizontal flows using particle tracking or PIV. 
Tubular bioreactor designs comprise arrangements of horizontal, vertical and/or inclined 
tubes in which cells are can be in transient turbulent or turbulent flow to maximize 
mixing and equalize light exposure (Grima et al. 2001, 2003; Garcia- Gonzalez et al 
2005; Chisti 2007). It is important to establish the most efficient flow rates for maximum 
biomass production with minimal energy consumption, and how the optical properties 
of bioconvecting suspensions (transmittance) change with the flow and couple to growth. 
Towards this end, insights into the transition to turbulence in the presence of cells will 
be invaluable, as will an understanding of how cells are arranged within a developed 
turbulent pipe flow (see Lewis 2003). The experiments presented here represent a first 
step in this direction. Theoretically, bioconvection in horizontal tubes has not yet been 
tackled; the main challenge is to predict the average inclination and speed of plumes, 
as well as the flow transitions we have observed, as functions of concentration and 
tube diameter. An open question is how to predict the effective transport properties of 
swimming cells in laminar and turbulent regimes in tubes of arbitrary orientation. 
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